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We study spectral properties of the Fokker-Planck operator 
that represents particles moving via a combination of diffusion 
and advection in a time-independent random velocity field, 
presenting in detail work outlined elsewhere [J. T. Chalker 
and Z. J. Wang, Phys. Rev. Lett. 79, 1797 (1997)]. We cal- 
culate analytically the ensemble-averaged one-particle Green 
function and the eigenvalue density for this Fokker-Planck 
operator, using a diagrammatic expansion developed for re- 
solvents of non-Hermitian random operators, together with a 
mean-field approximation (the self-consistent Born approxi- 
mation) which is well-controlled in the weak-disorder regime 
for dimension d > 2. The eigenvalue density in the complex 
plane is non-zero within a wedge that encloses the negative 
real axis. Particle motion is diffusive at long times, but for 
short times we find a novel time-dependence of the mean- 
square displacement, {r'^) ~ t^'"^ in dimension d > 2, associ- 
ated with the imaginary parts of eigenvalues. 

PACS numbers: 05.10Gg, 05.60.-k, 47.55. Mh, 02.70.Hm 



I. INTRODUCTION 

In this paper we study classical diffusion in the pres- 
ence of advection by a random, time-independent flow 
field, with an emphasis on spectral properties of the cor- 
responding Fokker-Planck operator. Our motivation is 
two-fold. First, the mathematical problem of calculat- 
ing the properties of a random, non-Hermitian differen- 
tial operator is interesting in its own right, and, as far 
as we know, has received little attention until recently. 
Second, spectral decomposition is a natural approach for 
investigating the diffusion-advection problem. Amongst 
other things, the spectrum contains information about 
the time-dependence of the effective diffusivity, and our 
results in fact reveal a short-time regime which appears 
not previously to have been discussed. 

Diffusion-advection problems can arise in a variety of 
physical settings, including turbulent diffusion of tracer 
particles in geological systems, the temperature field in 
Rayleigh-Benard convection, and flow through porous 
media; an extensive review is provided in the article by 
Isichenko |^. The interplay between advection and dif- 
fusion may either greatly enhance or strongly inhibit the 
long-time transport of particles, depending on the nature 
of the flow fleld. Broadly speaking, compressible flows 
inhibit transport, since fluid sinks at which flow lines 
converge will act as particle traps, while incompressible 



flows transport particles more effectively at long times 
than diffusion alone. In either case, the effective diffu- 
sivity provides a macroscopic measure of the combined 
consequences of molecular diffusion and of advection. 

Diffusion-advection problems have long been a subject 
of research in fluid dynamics. Early work dates back to 
that of Taylor |g] and of Richardson [^ , who proposed the 
notion of effective diffusivity. Subsequently, Batchelor 
analyzed one-point and two-point correlations in homo- 
geneous turbulence, relating the diffusion coefficient to 
the temporal correlation of the velocities Q. Recently, 
much analytic progress has been made in a special model 
considered by Kraichnan B , in which the velocity field is 
zero-mean Gaussian distributed and delta correlated in 
time, a special case for which the closure problem is ab- 
sent. In contrast, our interest in the following is in fiow 
fields that are time-independent. 

The effect of various forms of time-independent spatial 
disorder on classical diffusion has attracted considerable 
attention from statistical physicists, especially following 
the proof by Sinai that a type of random advection in one 
dimension results in dramatically sub-diffusive motion at 
long times B . Several groups of authors fTJ-p^ have in- 
vestigated scaling in these systems, particularly within 
the framework of renormalization group theory. Behav- 
ior is dependent on dimensionality, and long-time motion 
is diffusive only above an upper critical dimension, which 
is two if the advecting velocity field has only short-range 
correlations. While the focus of that work, reviewed in 
Refs. 11 IJ] has been long times and systems at or below 
the upper critical dimension, we restrict ourselves in the 
following to systems above the upper critical dimension, 
and consider transient as well as long-time behaviour. 

In addition to the study of diffusion-advection prob- 
lems as critical phenomena, many rigorous results, pro- 
viding bounds on the effective diffusivity and on the lim- 
iting behavior at long times, have been proved by math- 
ematicians using the methods of multiscale analysis and 
variational principles ||l^-[l5|] . 

Against the background of this varied literature, it is 
perhaps suprising that little seems to have been known 
until recently about spectral properties of the Fokker- 
Planck operator for random diffusion-advection prob- 
lems, even though the corresponding aspects of random 
Schrodinger operators have been studied very extensively. 
In fact, the spectrum of the Fokker-Planck operator de- 
pends very much on the nature of the velocity field re- 
sponsible for advection. In general, this velocity field 



can be separated into an incompressible (divergence- free) 
part, and a remainder, which can be expressed as the 
gradient of a scalar potential. Pure potential flow is a 
special case of some importance, especially since in one 
dimension any velocity field can be written in terms of a 
scalar potential. In this special case the Fokker-Planck 
operator is related by similarity transformation to a Her- 
mitian operator [|l6| , and therefore has a purely real spec- 
trum. Some previous work has built on this transforma- 
tion, showing that anomolous diffusion at long times in 
one dimension is connected with a logarithmic singularity 
of eigenvalue density of the Fokker-Planck operator [|7| , 
while other work has been concerned in two dimensions 
with the opposite limit of purely incompressible flow, an- 
alyzing the connection with the quantum random flux 
problem, and studying numerically the spatial decay of 
the Green function fig ]. 

More generally, as we show, the non-Hermitian na- 
ture of the Fokker-Planck operator is an obstacle that 
prevents one from transferring in a straightforward way 
the techniques developed for random Schrodinger oper- 
ators. For an unrestricted flow field, one expects eigen- 
values of the Fokker-Planck operator to occupy a finite 
area of the complex plane. As a result, the correspond- 
ing Green function is non-analytic throughout this area, 
and established perturbative approaches to calculating 
disorder-averaged Green functions, which depend on an- 
alytic continuation, are inapplicable. This difficulty has 
been faced previously in the study of certain ensembles of 
non-Hermitian random matrices, for which the eigenval- 
ues are known to be distributed uniformly within a circle 
in the complex plane |l^,Q. For these ensembles, spe- 
cial techniques have been developed |pl|,p2| , which do not 
generalise immediately to spatially extended problems 
such as the one we are concerned with. Recently, several 
groups independently p^-Eq] have found that one can 
make progress by constructing, from the non-Hermitian 
operator of interest, a Hermitian operator with a 2 x 2 
block structure, and applying standard methods to this 
enlarged Hermitian operator. We describe this approach 
below, and apply it to the Fokker-Planck operator, treat- 
ing the advection term as a perturbation to the diffusion 
term, a limit is known in the fluid dynamics as the small 
Peclet number regime. 

The general method described here might be applied to 
variety of problems involving random non-Hermitian op- 
erators in which behavior at weak disorder is of interest. 
A number of such problems have attracted recent atten- 
tion, including asymmetric neural networks [ pl[ and the 
statistical mechanics of flux lines in superconductors with 
columnar disorder p7| , and the Schrodinger equation for 
particles moving in a random imaginary scalar poten- 
tial [E8[. Recently, building on the formalism outlined 
in Sec. H, useful links have been established between ||29| ] 
problems of this kind and certain Hermitian localization 
problems. 

The present paper provides a detailed account of work 
presented previously in outline elsewhere p3| . In partic- 



ular, we describe for the first time calculations for a ran- 
dom flow field with arbitrary relative strength to the com- 
pressible and incompressible components. The remaining 
sections are organised as follows. In Sec. II, we define the 
problem and outline the Green function approach used 
in our calculations. In Sec. Ill we apply this approach to 
the Fokker-Planck operator for the diffusion-advection 
problem with a random velocity field. We obtain the 
Green function for the Fokker-Planck operator, and its 
eigenvalue density in the complex plane. We describe nu- 
merical tests of our results in Sec. IV, and present some 
technical details in two appendices. 



II. PRINCIPLES OF CALCULATION 

Particles that move by a combination of diffusion and 
advection have a density, n, which is a function of time, 
t, and position, r, and evolves according to the Fokker- 
Planck equation 



d^n = Cfp{n) = DV^n - V ■ (Yn) 



(1) 



where D is the molecular diffusivity, and V is the 
background velocity field, which we take to be time- 
independent and random. For definiteness, we apply pe- 
riodic boundary conditions to n(r,i) at the surface of 
a d-dimensional cube of volume O, but treat in detail 
only the limit SI — > oo. The Fokker-Planck operator, 
Cfp, is a sum of two contributions: the diffusion term, 
£0(71) — DS/^n, is Hermitian, while the advection term, 
Ci{n) — —V • (Vn), is not, and our attention through- 
out this paper is focused on their combined effects on the 
properties of Cfp as a random, non-Hermitian operator. 
Our choice for the probability distribution of the ve- 
locity field is intended to be the simplest: we take it to be 
Gaussian with zero mean, and with correlations of V(r) 
that are as short-range as possible. It is specified by the 
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where 



V(k) = (27r)-'' / d'^re 



-ikr 



V(x) 



(2) 



(3) 



(...) indicates an ensemble average, and Fi and r2 repre- 
sent the strengths of the incompressible and irrotational 
parts of the velocity field, respectively. For the diffusion- 
advection problem to be well-defined, it is necessary to 
impose a short-wavelength cut-off. A, on the spectrum 
of velocity fluctuations, so that Eq. (^ is supplanted by 
VQ(k) — for |k| > A. A dimensionless measure of the 
disorder strength involves a frequency scale, which we 
denote by uj: j,{uj) = {r,/\uj\'^){\u;\/Dy/^+^ is the ratio. 



raised to the power d + 2, of the distances travelled by a 
particle due to advection and to diffusion, in time u!~^. 
The fact that 7(0;) oc |w|(''^2)/2 identifies d = 2 as the 
upper critical dimension for the problem: in dimension 
(i>2, asa;^0,7(a;)^0 and one expects the long-time 
behaviour simply to be diffusive, while for d < 2, ^(lu) 
diverges as w — » 0, and long-time behaviour is dominated 
by advection. We restrict our discussion in this paper to 
the regime above the upper critical dimension. 

Our calculations center on the disorder-averaged one- 
particle Green function, which in the position-space and 
time domain is (n(r,t)), where n(r, i) satisfies Eq. (|l|) 
with initial condition n(r, 0) = 6{r). We obtain this from 
its Laplace transform. 
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(4) 



in which lu, in general complex, is the transform variable 
conjugate to t. It is useful to display the spectral de- 
composition of both these Green functions. Let {L\\ and 
\Rx) be the left and right eigenvectors of Cfp with eigen- 
value A. For finite system size, the spectrum is discrete, 
the eigenvalues are non-degenerate with probability one, 
and {(La|}, {|-Ra)} are complete, biorthogonal bases. In- 
troducing a plane-wave basis, (r|p) — Q'-^/^e^^''^, the av- 
erage spectral density is defined as 



C{p,u;)^{J2{p\R,){Lx\p)Sicj-X)) 



(5) 



Then 



(n(r, t)) = (27r)-'^ / d'^pe'^" / d^iue'^'C{p, w) . (6) 

The average Green functions are diagonal in the plane- 
wave basis, and we denote the diagonal elements oi g{uj) 
hy gP{uj): 



gn^) 



,.,Cip^ 

LU — X 



(7) 



where the dependence is only on p, the magnitude of p, 
in the large f2 limit, because of spatial isoptropy. Finally, 
the eigenvalue density is 



pM=J7"^^(5(cj-A). 



(8) 



Both C(p, A) and p(A) can be calculated from the Green 
function using the identity 



S{lu-X) = - ^ 



1 



■K duj* LO — X 



(9) 



with the notation d/duj* = {l/2){d/dx + id/dy) for lu = 
X + iy and x, y real. The eigenvalue density will turn out 
to be non-zero over a finite area of the complex uj plane, 
and not just on the real axis as would be the case if Cfp 
were Hermitian. 



Standard methods l^^l for calculating average Green 
functions for random Hermitian operators make exten- 
sive use of the fact that they are analytic in the com- 
plex cj-plane, except on the real axis. By contrast, g{Lo) 
is non-analytic throughout the (as yet unknown) area 
of the complex lo plane in which p{lu) is non-zero. Be- 
cause of this, a new approach is required. We first embed 
the original non-Hermitian operator, in the combination 
LU ~ C = A, in a Hermitian operator Ti which has twice 
the dimension of £, setting 



Ti — Hq + Til 



(10) 



where |^1| 
Hq = 



77 
-r/ 



and Til = 



A^ 



A 




(11) 



The inverse of H exists for real 77 ^ 0. Its average is 



G 



Gil 
G21 



G12 
G22 



_ ( (77[r?2+yl^t]-i) 
" {{A^lri' + AA^-^) 



{A[v^ + A^A]-^) , 

^2 ' A^A]-^) '■ ^^^^ 



{-vVi^ 



Note that the desired Green function is obtained using 
g{iju) = \va\ri^QG2i- Since Ti. is Hermitian, we can cal- 
culate G, and hence (/(w), with extablished techniques. 
Specifically, G can be expanded as a power series in Tii 
and in Go = Ti-Q^ 

G = Go — {GqTLiGq) + {GqTLiGqTLiGq) + .... (13) 

which leads to a Dyson equation with self-energy E 



G — Go + Go^G . 



(14) 



Working in the plane- wave basis, G and E each consist 
of four diagonal blocks, with elements Gij{p) and S^ (p), 
for i, j — 1 OT 2. It is convenient to represent the power 
series diagrammatically. As shown in Fig. (nf), we have 

two propagators, G\i = l/i] and G22 = —1/??, denoted 
by a single and double line, respectively, and two vertices, 
A and A''. It is natural to separate each vertex into a 
constant part and a random part: case oi A, cu — Co and 
~Ci. LU ~ Co is diagonal, with diagonal elements LUp = 
LU + Dp^ , and averages of the random parts are specified 

by the four cumulants: (CiCi), (CiC\), (C[Ci), and 

(44). 

The leading contributions to G at weak disorder are 
resummed by the self-consistent Born approximation 
(SCBA). This approximation is defined by the expression 
for the self-energy shown diagrammatically in Fig. (0), in 
which the internal propagators represent the full Green 
function G. We demonstrate in appendix A that cor- 
rections to the SCBA self-energy are small in powers of 
7(1^1). 
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FIG. 1. Propagators and vertices. 
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FIG. 2. Diagrammatic representation of disorder average. 
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FIG. 3. Full Green function and self-energy diagrams in 
the SGBA. 

There are useful interrelations among the four blocks 
of Gij , and also among those of S^j . Quite generally, in 

the limit e ^ 0, one has G21 = GJj and S21 = '^W: us- 
ing the fact that each is diagonal in the plane- wave basis, 
one also has T,ii{p) = — E22(p) and Gii(p) = — G22(p). 
It is therefore sufficient to calculate only Gii(p), G2i(p), 
Si2(p) and S22(p)- They are related via the Dyson equa- 
tion, which yields 



Gn{p) = [v + ^22{p)]C{p) 



(15) 



G2i{p) = [-I:Up)]C{p) 



where 



c{p) = [\j:i2{pr + \v+^22{p)n 



2i-l 



(16) 



(17) 



In turn, the SCBA generates expressions for Y.ij{p) in 
terms of Gij (p) : referring to Fig. (^ , we have 



Si2(p) 



d''q{Cl{p,q)Cl{q,p))G2i{q) (18) 



^22{p)=J dSl^C\{p,q)C,{q,p))G^M- (19) 

The solution to this set of equations, with the velocity 
correlations of Eq. (||) , simplifies considerably in the spe- 
cial case Fi = r2, for which the calculations have been 
described previously |2^. The solution in the general 
case is presented in the following section. 



III. APPLICATION OF SCBA 

The expressions for the self energy within the SCBA, 
Eqs. (|l^) and ( [l9| ) become, with the Gaussian velocity 
distribution of Eq. (0) , 



.A 

Si2b) - -cc-p + Fi / dSiii ■ q)I]*2(g)G(q) + 
Jo 

^l2iq)C{q) 



(r2_ri)/Al^(^-P)][P-(^-P 



|q-p| 



(20) 



.A 

S22(p) = El / d'^qq^irj + S22(Q)]G(g) + 

JO 

(T^-r,) f dU^^jP^Z^[V + ^22{q)]C{q) (21) 



|q-p| 

where the fact that /^ d'^q{p.q)'E,l2iQ)(^{'l) ~ fo'' sym- 
metry reasons is an immediate simplification. Eqs. (Efl) 
and (pi|), with Eq. (|3); are coupled integral equations, 
and presumably in general not exactly solvable. To make 
progress, we treat weak disorder, for which the problem 
can be reduced to a system of algebraic equations. As 
a starting point, consider the form of the function C{q), 
Eq. (p^), which appears in the integrands of Eqs. ( pO[ ) and 
(pi|). First, in the absence of disorder, 'Ei2{q) = —iOq, 
S22 — and for 77 ^ 



C{q)^ 



1 



{x + Dq^Y + y2 



(22) 



and so, if oj = a; -f- iy lies on the negative real axis, C{q) 
has a divergence as a function of g, at q — qo, where 
x + DqQ = 0. In the presence of weak disorder, we antici- 
pate that G{q) will be large only for x < 0, \y\ sufficiently 



small and q ^ qq. We also expect that, if to lies in the 
right half of the complex plane, or if \y\ is large, so that 
C{q) is small throughout the range of integration over g, 
then a good approximation for Sjj(p) should be obtained 
from evaluating Eqs. (EQ) and (21) by iteration, using the 
disorder-free form for C(g) in the integrands. This ap- 
proach is analogous to the first-order Born approxima- 
tion of scattering theory. Conversely, if a; < and \y\ is 
small, we expect that it will be necessary to determine 
C{q) self-consistently, but with the simplification at weak 
disorder that the dominant contribution to the integrals 
of Eqs. (^0|) and ( [21[ ) is from the shell of wavevectors on 
which g « go- For a; < 0, we therefore evaluate the 
integrals by replacing C{q) with I{uj)5{q^ — g^), where 
/(w) = /p d{q''^)C{q'). This procedure yields the leading 
behaviour at small Fi, r2, provided the integrals involved 
are convergent at small and large g, which is the case in 
dimensions 2 < d < 4. With this approach we find 



Si2b) = 



(r2 



ri)^s 



[n ■ {n -rfo)]Ko • (r 



i2(go)/(w) X 
"o)] , 



S22b) = ^[r7 



S22(go)]/(w) X 



-)n 



il + (i2 -ilj( — 77 rr^ — )fi 



(23) 



(24) 



where Sd is the surface area of a unit sphere in d dimen- 
sions, n and no are d-dimensional unit vectors in direc- 
tions corresponding to those of q and p in Eqs. ( po| ) and 
(pi|), and {. . .)n denotes an angular average on fi. These 
averages have the values 



[n ■ (n - no)] [no • (n - n'o)] 



(n - no)^ 

[n • (no 



(- 



h)Y 



yn-no) 



/" — r, 



(25) 
(26) 



Notation is simplified by introducing the new variables: 

r; = T,Sdqt/2, A = (r; - r;)/2 and r = (rU t'^)/2. 

In these terms, the solution to Eqs. ( p3| ) and ( |24| ) is 



Si2(p) = - [x + Dp^ +i 
^^' 1-TI{lu) 



y 



1 - A/(cj) 



(27) 
(28) 



Substituting these expressions into Eq. (|l7|), we obtain 

C{q) = [{x + Dq^f + l3^]-\ (29) 

where 



/3^ 



1 - A/(cj) 



^1 



1 - TI{uj) 



(30) 



and hence (for x ^ —(3) 



I{lo) 



diq')C{q) 



d{q')C{q)^—. (31) 



Thus, finally, /(w) satisfies 



yi{^) 



1 - A/(w) 



ril{uj) 



l-TI{uj)J \D 



r IT \'^ 



(32) 



The behavior of the solution to this equation in the 
limit 77 ^ depends on the value of y. If y is sufficiently 
large {y > ys, where ys is determined below), r/(w) < 1 
and the second term on the left of Eq.,(p2|) vanishes as 
7/ — > 0. The limiting solution is then simply 



Iiu;) = 



Dy + TrA 



(33) 



In that case 



gP{Lu) ^ l[inG2i ^ \Dp^ +x + i{y + An/ D)] \ (34) 

(35) 

This expression is manifestly analytic in uj, implying from 
Eq. (g) that p(cj), the eigenvalue density of the Fokker- 
Planck operator, is zero in this part of the w-plane. A 
similar analysis applies for y < —ys- In contrast, for 
\y\ < yB, the second term on the left of Eq. (p2) remains 
non-zero as 77 — » and the limiting solution is 



lico) = r 



(36) 



The boundary between these regimes is evidently \y\ = 
yB with 



VB 



= ^ = ^7i(-)-Nocri|.r. 



D 2 

Inside this boundary, 

Dp2 + 2- _ iyT/r[ 



/(^) = 



{Dp^ + xY + [ttT/DY 



(37) 



(38) 



which, in the same spirit as our representation of C{p), 
may be written as 



5» = -^5(p^-gg). 

J- 1 

From this we obtain, for a; < and \y\ < yB 

D 



p{uj) 



(27r)'^+iFi|5 



(39) 



(40) 



Finally we turn to the time-evolution of the particle 
density. From g^icu) we calculate 



C{p,^) = ^6{p'~ql) 



(41) 



Using Eq. (0), we find 



(n(r,t)) 



1 



(27r)'i 



d^pe 



jp.r 






(42) 



where loq — ys, evaluated at x = —Dp^: ujq = 
7TSdTip'^/2D. The integrand has two factors: the first 
one, e^^P *, would arise simply from diffusion; the sec- 
ond, sm{LL)Qt)/uJot, is a consequence of advection. At weak 
disorder, when the SCBA is a controlled approximation, 
the advective factor differs significantly from 1 only where 
the diffusive factor is small, so that (r^) ^ Dt, indicat- 
ing the normal diffusive behavior. By contrast, at strong 
disorder (when the SCBA is simply a mean-field approx- 
imation), it is the advective factor that sets the width of 
the density profile at short times, and (r^) ~ (Tit/D)'^^'^. 



IV. COMPARISON WITH NUMERICAL 
SIMULATIONS 

As a test of our analytical approach, we compare re- 
sults for the boundary of the support of the eigenvalue 
density with numerical calculations of the eigenvalue dis- 
tribution. The numerical calculations involve diagonalis- 
ing a discretised version of the Fokker-Planck operator, 
as set out in this section. In order that the test of an- 
alytical reults is as direct as possible, we adapt the an- 
alytical theory for the discretised operator, in the way 
summarised in Appendix B. We can make the compar- 
ison for a two-dimensional system despite the fact that 
the theory developed in Sec. Ill has long-wavelength di- 
vergences in dimensions d < 2, since the system sizes we 
study are small enough to cut off the divergences. 

We discretise the Fokker-Planck operator keeping in 
mind the conservative form of the Fokker-Planck equa- 
tion 



dn 
J = DVn - Vn . 



(43) 
(44) 



On a two-dimensional, square lattice with coordinates 
i,j, velocity components Ui^ and Vi,j are defined on hor- 
izontal and vertical links, respectively, and densities riij 
at nodes. Using center difference, the discrete eigenvalue 
equation becomes 

Xriij = [D- -u^+ijrii+i^j + [D + -u.^_ij]ni-i.j 



-AD 



-(wi+i J - u,,_i J + w;.j+i - u,j_ 1 )]n,j (45) 



We restrict ourselves to the case Fi = r2 = F, for which 
the velocities (it, v) are Gaussian random variables with 
zero mean and variance (27r)^F. A cutoff |m| < 2D is im- 
posed to ensure that matrix elements of the discretized 



operator are non-negative, a requirement for the transi- 
tion matrix in any Markov process, which in turn implies 
that the real parts of the eigenvalues are necessarily non- 
positive. The eigenvalues found by numerical diagnaliza- 
tion, using 50 realisations of a 32 x 32 lattice with D = 1 
and (u2) = (t,2) ^ 0.25 so that F = 0.25(27r)-^ are 
shown in Fig. (Q), as well as the boundary to the eigen- 
value density calculated for the discretized problem on 
a lattice of the same size (full line), and for the contin- 
uum theory (dotted line). The full line is obtained, fol- 



lowing Appendix B, as the lo values which solve Eq. B8 . 
The dashed line is calculated from Eq. B9: it has slope 
Tt:S2/2 = 1/16. 
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FIG. 4. Distribution of eigenvalues in the complex plane 
(dots) and calculated boundary to p(a;) for the lattice theory 
(full line) and continuum theory (dotted line) . 

The finite lattice spacing manifests itself in a lower 
bound to the real part of eigenvalues, and the finite sys- 
tem size studied results in oscillations in the position of 
the boundary that are apparent in both the numerical 
and analytical results. A further finite-size effect, clear 
in the data but not captured in the SCBA we have pre- 
sented, is that a finite fraction of eigenvalues are purely 
real: phenomena of this kind in random matrix problems 
have been analysed in Refs. ( pj) and ( Urn). 



V. SUMMARY 

We have described in detail a general technique for 
calculating average spectral properties of non-Hermitian 
random operators. We have applied this technique to 
the Fokker-Planck operator that represents particle dif- 
fusion in the presence of random advection. In this way 
we obtain for the first time the eigenvalue distribution 
of this operator in the complex plane. The approach is 



formally exact in the weak disorder limit, and numer- 
ical tests show that it is remarkably accurate even for 
moderately strong disorder. Finally, we have inferred 
the time-dependent effective diffusivity from the ensem- 
ble averaged Green's function, obtaining a novel scaling 
in its transient behaviour in the strong disorder limit. 
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APPENDIX A: CORRECTIONS TO SCBA 

We demonstrate in this appendix that corrections to 
the SCBA expression for the self-energy are small in the 
coupling constant, 7(|aj|), considering for simplicity the 
case Fi = r2 = F. We compare the SCBA self-energy, 
5^11 (p), with the contribution arising from the diagram 
illustrated in Fig. (^, which we denote by T,{i{p). 




V(2) 



FIG. 5. Leading correction to the SCBA self-energy. 

(2) 

Our goal is to compute the ratio T,\-^ {p)/T,ii{p) in the 
weak disorder limit. We begin by calculating the correc- 
tion 

sf^p) = r2 / / AiA2(p-q2)'x 

Jq Jo 
G22(gi)Gii(q2)G22(p + q2 -qi). (Al) 

First, consider the region of the complex to plane in 
which the Fokker-Planck eigenvalue density, calculated 
within the SCBA, is zero (all y for a; > 0, and \y\ > ys 
for X < 0). In that region, for small rj, Gulp) oc ry, 

I]ii(p) oc T], 'Sl^ (p) oc rf and the ratio has the limiting 
value T,\^'{p)/T,ii{p) = for 77 ^ 0. Second, consider 
the complementary region (a; < and \y\ < ys)- In thi s 



region, we use the same approach for evaluating Eq. ( Al ) 
as was applied to Eqs. (23) and (pj) in Sec. III. Substi- 
tuting for Gii{p) in terms of C'(p)7^replacing C{p) with 
I(u))S{p'^ - Qq) where Dqq -|- x = 0, and using Eqs. ^^) 
and ( p2[ ) to obtain Sii(p) and I{oj), we find, for p = qa 
in the limit rj —> 



(2)^^^ _ (vb - y') 



^YAp) - 



2 _,,2^3/2 



r<z§ 



-/ 



(A2) 



where / is an angular average over directions of qi , q2 , 
represented by the unit vectors rii and 712 

/ = {{n\ ■ n2)^S{[n + n^ - ni]^ - l)),rf,,f, = 0(1) . (A3) 

Thus, for \y\ < yB, corrections to the SCBA have relative 
size 



and, as advertised, are small at weak disorder. 



(A4) 



APPENDIX B: STUDY OF DISCRETIZED 
PROBLEM 

On a lattice, with site labels a, /? in arbitrary dimen- 
sion, the advection part of the Fokker-Planck operator, 
—V • V, is represented by a matrix with elements 



A. 



1 



Q/3 



:Ua[3 



(Bl) 



for a, (3 nearest neighbor sites. It satisfies Aap — —Af^a 
and has diagonal elements 



A, 



^^oH""^ 



(B2) 



so that X^a^C"^'/^) ~ ^' ^^ required for conservation of 
particle density. Other elements are zero. We take Aap 
to be Gaussian distributed with zero mean, and variance 

{AanAap) - rvr^ (A^pAp^) = -Ttt^ (B3) 

(^aa^a/3) = Ttt^ (A^a^Ap^) = -Ttt^ (B4) 

(^L) = 4F7r2 (A^^App) = ~rn\ (B5) 

for Q, f3 nearest neighbor sites. With the plane-wave 
basis (rajk) — jv^^/^e**''"", where Tq, is the position co- 
ordinate of site a on a lattice with N sites in total, we 
need ((q|yl|k)(k|A-'"|q)), which takes the value 

{{q\A\k){k\A^\ci)) = 

ST^ I A A \ iq-(ri-r3)+ik-(r4-r2) 

7VT2 / J \^rir2^r3r4/ '^ 



7V2 



rir2r3r4 

4^2r 



N 



[(1 -cosga;)(l + coskx) + 



(1 — cosgy)(l -I- cosky) 

Similarly, 

((q|^|k) (k|^|q)) = 
(27r)2F 



(B6) 



N 



sin(fc^)sin((7j;) +sin(fc;,)sin((jj^)] . (B7) 



Repeating the style of calculation described in Sec. Ill, 
we find that the support of the eigenvalue density of the 
discretized Fokker-Planck operator has a boundary de- 
termined by 



(27r)^r Y^ {1 + cosqx){2 — cosqx — cosqy) 



N 



E 



\UJ - 



2D[2 — cos qx — cos qy] 



= 1. 



In the continuum limit we recover the result 

„2 



A 

l^^- 



Dq 



.212 



1 



(B8) 



(B9) 



and hence Eq. (t5 
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